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Abstract 



Semiclassical trace formulas are examined using phase space path 
integrals. Our main concern in this paper is the Maslov index of the 
periodic orbit, which seems not fully understood in previous works. We 
show that the calculation of the Maslov index is reduced to a classifi- 
cation of connections on a vector bundle over S 1 with structure group 
Sp(2n, R). We derive a formula for the index of the n- repetition, and 
show that a Bohr-Sommerfcld type quantization condition including 
quadratic fluctuation around the orbit is derived using this formula. 
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C Useful formulas 



1 Introduction 



The semiclassical trace formula developed by Gutzwiller Q is one of 
the most important tools to study the spectrum of non-integrable sys- 
tems. However, the original derivation of this formula is very com- 
plicated, and the canonical invariance of the formula (especially the 
Maslov index) is unclear. 

In this paper, we examine the semiclassical trace formula using 
phase space path integral to clarify the canonical structure of this 
formula. Our main concern here is the Maslov index, which is an 
additional phase factor appearing in the trace formula. 

Although many studies have been made on the geometrical prop- 
erties of this index, it doesn't seem to be thoroughly understood. For 
hyperbolic orbits, Robbins B showed that this index is equal to twice 
the winding number which is defined by the invariant manifolds around 
the orbit. In this case, the Maslov index of the n-repetition of the orbit 
fx n is equal to n/xi. However, Brack and Jain J3| investigated periodic 
orbits in anisotropic harmonic oscillator (these orbits are elliptic), and 
showed that [i n is not equal to n\i\. This result suggests that we need 
new geometrical picture to understand this index. 

We start with the phase space path integral of the partition func- 
tion: 

Z(T) = JvpVqexp ^ I (pdq - Hdt) . (1) 

The density of states is obtained as the Fourier-Laplace transformation 
of the partition function: 

p(E) = --Img(E + ie), (2) 



9(E) = (3) 
h/ — H 

= _ / dTe lET/h Z(T). (4) 
in Jo 

We obtain the semiclassical trace formula by applying the stationary 
phase approximation to the path integral. 
The stationary condition reads 

5 <j)(pdq - Hdt) = 0, (5) 

which leads to the Hamiltonian equation of motion. Since the periodic 
boundary condition of the path integral (Q) has the effect of closing 
the paths, the solutions of (pi) are periodic orbits. 
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Thus the partition function is approximated by a sum over the 
periodic orbits: 

(6) 



Z(T) = ^^exp 



p.O 



Here, R = § pdq — Hdt is the classical action of the periodic orbit, and 
K is the contribution of the quadratic fluctuation around the orbit: 

K = J Vxexp[i5 2 R[x(t)\], (7) 

x(t) = -j=(6q,6q). (8) 

The Maslov index appears as the phase factor of this quadratic path 
integral. Since this is a kind of Fresnel integral, the phase is determined 
by the signs of diagonal elements of the quadratic form S 2 R. 

Strictly speaking, we should write (|l|) in a discrete form to obtain 
well-defined continuum limit. Since the way of the discretization re- 
flects the operator ordering of the Hamiltonian, we should treat this 
problem carefully Q. In this paper, we always use the mid-point pre- 
scription, which corresponds to the Weyl-ordering of the operators. We 
can formulate the canonical transformations of the path integral most 
clearly in this prescription, as we will see in §|3| 

The central idea of this paper is that the Maslov index of the peri- 
odic orbit is determined by the linearized symplectic flow around the 
orbit. The set of displacement vectors {x(t)} is considered to be a vec- 
tor bundle over S 1 , and the flow around the orbit define a connection 
(Fig. [I]). Mathematically speaking, our problem is the classification of 
connections on the fibre bundle. However, the quadratic path integral 
around the periodic orbit is not invariant to all canonical transforma- 
tions. If the canonical transformation is topologically non-trivial, the 
path integral may change the sign. This is essentially the same as 
the global anomalies of the gauge field theories Jl(], [□]]. Therefore we 
regard two connections as equivalent if they are connected by a topo- 
logically trivial canonical transformation, and classify the connections 
by this equivalence relation. 

This paper is organized as follows. In we review the derivation 
of the phase space path integral of the partition function and the sta- 
tionary phase approximation. Note that the discussion in this section 
is not restricted to the time-independent Hamiltonians. §^ is the main 
part of this paper, and we treat the quadratic path integrals with peri- 
odic boundary conditions generally. We calculate the path integrals by 
reducing them to simple normal forms using canonical (gauge) trans- 
formations. We also derive a formula for the Maslov index of the 
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periodic orbit 



Figure 1: The set of displacement vectors is regarded as a fibre bundle, and 
Hamiltonian symplectic flow defines the connection. The structure group of 
this space is Sp(2n, R) 



n-repetition of the orbit. In §£|, we discuss the trace of the resolvent 
g{E) for time-independent systems. If the Hamiltonian of the system is 
time-independent, the periodic orbit in it have at least one zero-mode. 
Therefore we first discuss the integration with respect to zero-modes. 
Then we execute the Fourier transformation of the partition function 
and derive the semiclassical approximation to g{E). We also show that 
the sum over repetitions of a orbit leads to a Bohr-Sommerfeld type 
quantization condition. Similar attempts have been made by many 
authors, but our result is the first one which includes the effect of the 
Maslov index correctly. 
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2 Path integral and stationary phase ap- 
proximation 

2.1 Path integral representation of the partition 
function 

In this subsection, we derive phase space path integral of the partition 
function 

Z(T) = TW(T). (9) 
Here, U (T) is the time evolution operator 



U(T) = Pexp 



i 

h Jo 



T 

dtH(t) 



(10) 



and H(t) is the Hamiltonian of the system. P denotes the path ordered 
(or time ordered) product. 

U(T) can be decomposed into a product of infinitesimal time evo- 
lution operators: 

U(T) = f[{l- UtH( tj )\+0(l/N), (11) 
j= i L > 

where At denotes T/iV, and tj denotes jAt. Therefore Z(T) can be 
written as 

Z(T) = Jdq(q\U(T)\q), (12) 

- N 

= ^jU^jW-jrteHiti)^). (13) 
i=i 

Here, the variables {qj} are cyclic and qi is considered to be the same 
as qj if i = j(modN). This notation is used throughout this paper. 

The Feynmann kernel for the infinitesimal time interval can be 
written as 

<«ii - iA<ij( ( ,)| 5j _,> . / ^ {i - iMH W („„ iiip+^M^W) 



J T^kt exp n [ Pj{qj Qj - l] Hw [ p » — 2— ' 'V 4i 



where n is the number of degrees of freedom, and Hw is Weyl transform 
of H: 

H W (p,q) = Jdve^ v l h {q-^\H\q+^). (16) 
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Therefore we obtain the phase space path integral representation of 
the partition function: 



A' 



Z(T)= hm TT 



3=1 



exp 



N 



Y, {P, ,( lj - - &tH(Pj> 

3=1 

(17) 

Here, we denote Hw as simply 7J, and q^ denotes (q^ +qr,-_ 1 )/2 (mid- 
point prescription). 



In the continuum limit, we denote this path integral as 
Z(T) = J VqVpexp i ^"(pdq - i?dt) 



(18) 



Note that the c-number Hamiltonian is defined by the Weyl transform 
of the q-number Hamiltonian in this paper. This c-number Hamil- 
tonian can also be obtained by arranging the q-number Hamiltonian 
into Weyl ordering, and substituting c-number variables for q-number 
variables Q. 

If the Hamiltonian is given as 



(19) 



integration with respect to p in ( |17j ) can be done easily, and we obtain 
the path integral representation of Z in coordinate space: 



z <n=x* n/( 



3=1 



/ m 



,2nihAt 



n/2 



dq, cxp 



. N 

3=1 



m ( q, -q, i 



At 



-V(q r t,) 



In the continuum limit, this formula can be denoted as 



Z(T) 



Vq exp 



- $ Ldt 



(20) 
(21) 



where L is the Lagrangian of the system. 



2.2 Stationary phase approximation 

Let us evaluate the path integral ( |l7| ) by the stationary phase approx- 
imation. We define a discretized action Rjy as 

N 

Rn(Pi,P 2 , -iPiv.tfn -iIn) = {Pj^j ~ «3'-i) ~ AtH (P 3 ^Qj^j)} ■ 

3=1 

(22) 



9 



Then the stationary phase condition reads 



= 



OR 



■N 



dp, 



= q ] -q^ l -^t—{p v q v t 3 ), 



(23) 



„ dR N At (OH , . dH. A. 



which leads to the Hamiltonian equation of motion in the limit N — > oo: 

(26) 



Since variables in iijv are cyclic, the solutions of the stationary phase 
condition are classical periodic orbits. 
The second derivatives of Rn are 



d 2 R 



N 



dpidpj 

■N 



a 2 R 



(27) 
(28) 

dPidqj ~ l " M;!JJ - (29) 

Here, H pp , H pq , H qq and A are Nn x Nn matrices, which are defined 



dq l dq J 



At{Hpp)ij : 

= -At(H qq ) itj , 

= &t{(A)ij — (H pq ) it j}. 



as 



(Hpp)i,j — 
{Hpq)i,j — 
{Hqq)i,j — 



Therefore 



(Pj,Qj,tj)Si 



d 2 H 

1 f d 2 ff 

2 \ <9p9q 

32 



(Pj+l)9j+l)*j+l)^i,j+l + 



d 2 tf 



1 f <9 2 iJ 

4 ("^2"^'+ 1 '^'+ 1 '^+ 1 )^'+ 1 + ^ 
- Sij +1 )/At. 



(30) 



(32) 

(33) 
(34) 



Z(T) = ^exp 
= H ex P 



-it 



p.o. 



■cl 



-rRcl 

h 



i^ooi (2^ 6XP 

e »f (M+,JV-M-,Jv) 



2ft 



lim 



^Idct^ 



(36) 



N 
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where 



<j>(PciQci - H(p cl ,q cl ,t))dt, 
(Spi,8p 2 , -,Sp N , Sq 1 ,Sq 2 , ....,Sq N ), 



At 



ffpp 

A T - HL 



A -H p 

—Hnn 



(37) 
(38) 
(39) 
(40) 



and n+ t N (a*-, at) is the number of positive (negative) integers. Since 
A becomes differential operator in the limit N — > oo, Z can be written 
as 



Z(T) = 



Tt. ^Det^ 



: exp 



-Rc 



Here, the differential operator Tj) is defined as 

P =Jj t -H"(p cl (t),q cl (t),t). 
J and H are 2n x 2n matrices defined as [] 

J 



H 







(-° 


o)' 


d 2 H 


d 2 H 


dp 7 

a 2 H 


dpaq 

a 2 H 


dqdp 


dq 2 



(41) 



(42) 



(43) 



(44) 



The functional determinant of If) is defined as the limit of the de- 
terminant oijp N : 

Det$> = lim dct^V (45) 

N^oo 

/i is the Maslov index of the periodic orbit, which is defined as 



lim 

N— >oo 



(46) 



Since If) N is even dimensional matrix, — (J-+.N is also even and fj, 
is integer if there is no zero-mode. 



1 Note that the definitions of Ip and J are different from those in the reference Q. (The 
signs are opposite.) Hence the definition of the Maslov index ( fl6| ) is also different from 
that in We follow the definitions in || in this paper. 
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2.3 Relation between Maslov index and Morse in- 
dex 

Let us consider the case in which the Hamiltonian of the system is 
given by 

H=^- + V(q,t). (47) 

In this case, the second variation of action around the periodic orbit is 

1 N ( At 1 At N 

-X T N ip N X N =-Y,\ 5 PA 5 9j \-- E VijQiSqj, 

(48) 



,2m r] ^ -*J-±<\ 2 

j=l v J i,j=l 



where V, „■ is dehned as 



a 1 \ d V d V 

V li = 4{-^(5j+i.*j+i)*ij+i + -^2(5i+i.*i+i)*iJ ( 49 ) 



<9 2 y _ d 2 v _ l 

+ 9 q 2^r t j) 6 i,j + -ggTilpW^j-i >■ (50) 
This quadratic form can be transformed as 
1 N At 2 

= -Es( ip i-S ( ^-M (51) 



2 V At / 2 



N 

= - E + T E I ) <I- ^i J<1. 6 <Ij- (53) 

Here, = Jp^ — m(<5q J - — <5q J _ 1 )/At, and the second term of RHS is 
the second variation of action in coordinate space path integral (po|). 
In the limit iV — > oo, Dq n reduces to a differential operator: 

D q = -m^ + v"(q cl (t),t). (54) 

(Here we used the partial integration (dq/dt) 2 — > —d 2 q/dt 2 .) We de- 
note the number of negative eigenvalues of iVn x iVn quadratic form 
Dq n as /iq,-,Ar. Then the Maslov index of the orbit is 

V = J lm o ' ( 55 ) 

iv— >oo Z 

= lim Nn ~« Nn -HQ.-*} (56) 

= lim Hq-.N- (57) 

N—>oc 
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phase space 


coordinate space 


partition function 


/ VpVq cxp j- §pdq - Hdt 


/ Vq cxp f § Ldt 


second variation 


V=J±-H" 


D q = -m^ + V" 


Maslov index 


(fi- - fi+)/2 


H- (Morse index) 



Table 1: phase space v.s. coordinate space 



This is the number of negative eigenvalues of Dq (Morse index) , which 
is finite if m is positive. Therefore the Maslov index is the Morse 
index in coordinate space path integral. We summarize the correspon- 
dence between phase space path integral and coordinate path integral 
in Table |. 
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3 General treatment of quadratic path in- 
tegrals 



In section |2|, we evaluate path integral by the stationary phase approx- 
imation, and it is necessary to calculate quadratic path integral like 
( |35| ) . In this section, we study the general structure of the quadratic 
path integrals. First we note that quadratic path integrals are inde- 
pendent of h. (We can show this explicitly by changing the variable 
Xn — > VTIXn in (|35|).) The period T can also be normalized by 
t —> Tt. Therefore it is enough to consider the case in which h = T = 1. 



3.1 Quadratic path integrals 

Let us consider the following quadratic Hamiltonian: 

1 T 

H{p,q 7 t) = -x h(t)x. 



(58) 



Here, h(t) is a symmetric 2n x 2n matrix, and x = (p,q). Since this 
Hamiltonian is already arranged into Weyl ordering, the corresponding 
classical Hamiltonian is 



H(p,q,t) = ^x T h(t)x, 



(59) 



where a; is a classical 2n dimensional vector (jp, q) . Our concern is to 
calculate partition function 



Z = TrU(T = 1), 

l-T 



U{T) = Pexp 



H(t)dt 



Z can be represented by path integral as in 



Z = lim 



lim 



(2n) Nn 



exp 



-X T N Ip N X N 



where 



n^oo ^det^l' 



N 



-H, 



pp 



N \ A r -H, 



p'i 



A - H, 
-H„ 



p'i 



(60) 
(61) 

(62) 
(63) 

(64) 
(65) 
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(Hqq)i,j,k,l 



— {h(tj + i)k,n+lSij + i + h(tj)k,n+lSij} , 



{h(tj + i) n +k,n+l{^i,j + l + Sij) + h(tj) n+ k,n+l(Sij + Si, 



A 



i,j,k,l 



(1 < i,j < N,l < k,l < n). 



(70) 



The quadratic path integral (|35|) appearing in the stationary phase 
approximation is the special case of (163) in which 



(71) 



h(t) = TH (p cl (tT),q cl (tT)). 
In the continuum limit, we write this path integral as 



Z = 



Vx cxp 



-x T ]px 



P = Jj t -h{t), 
x = (p,q). 



(72) 
(73) 

(74) 
(75) 



Our purpose in the following subsections is to calculate ( |73|). However, 
It is difficult to calculate (f73j) directly from the definition (|63|) . We will 
explain our strategy to calculate (73) in the next subsection. 



3.2 Outline of the discussion 

Our strategy to calculate (|73| ) is to transform the Hamiltonian matrix 
h(t) into a simple normal form by canonical transformations. 

Let us consider a time dependent symplectic matrix S(t). This 
matrix satisfies 

S T (t)JS(t) = J, (76) 
and we assume that S(t) satisfies periodic boundary condition: 

5(0) = 5(1). (77) 

S(t) defines a canonical transformation 

x\t) = S{t)x{t). (78) 
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The Hamiltonian matrix h(t) is transformed by S(t) as 



h^iS-'fhS- 1 - JS^- (79) 

We can see this transformation from a different point of view. As 
it was shown in ||, the set {a;(£)|0 < t < l,x(0) — x(l)} is considered 
to be a vector bundle over S 1 , and A(t) = Jh(t) is a connection on 
this bundle. Therefore canonical transformation defined by S(t) is 
considered to be a gauge transformation. Actually it is easy to verify 
that A = Jh is transformed by 5 as 

A — > SAS^ 1 + S-- — . (80) 
at 

The operator Ip can be rewritten as 

]p{t) = JD{t), (81) 

D(t) = j t +A(t). (82) 

D is a covariant derivative, and parallel transport defined by this co- 
variant derivative is the same as the Hamiltonian equation of motion: 

- + A(t)| x(t) = — (i) + Jh(t)x(t) = 0. (83) 

We can transform given A(t) into simpler form by gauge transforma- 
tions. However, the path integral ( |73|) is not invariable to all gauge 
transformations. 

Let us explain this point by a simple example, namely the time 
independent one dimensional harmonic oscillator: 



H=-(p 2 +q 2 ). (84) 



In this case, 



and 



a 
a 



—a d/dt 
—d/dt —a 



(85) 
(86) 



Eigenvalue equation 
is easily solved as 



)x n (t)=e n x n (t) (87) 
i n = -a + 2nn, 
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Figure 2: Flow of eigenvalues of Ip . 



, , / cos 2mrt W - sin 2rari \ . . 

Bb W = ^ s in2mrf )'\ cos2mrf J' (89) 

Note that each solution is doubly degenerate. If a = 0, H = and 
the Maslov index is considered to be in this case. As a changes 
continuously, eigenvalues {e n } also changes continuously and two of 
them change their sign when a crosses multiples of 2tt (Fig.|f). Thus 
we obtain the formula for the Maslov index: 

[ 2m (a = 2m7r) 

(See also appendix |b|). 

The canonical transformation defined by 

q /.\ _ f cos2fc7rf — sin2fc"7ri \ . . 

bk{ > = \ sin2/brf cos2fc7r< ) [J ) 

changes the Hamiltonian as 

fb 2 + 9 2 )-^^V + g 2 ). (92) 

Therefore the Maslov index changes as 

H^fi + 2k, (93) 

though the spectrum of Ip is unchanged. 

This relation means that the sign of the path integral changes if k 
is odd. Let us discuss this result from a different point of view. 

The classical time evolution operator of this system is represented 
by the matrix 

/cos at -sinai \ 
V sin at cos at / 
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On the other hand, the quantum time evolution operator U (t) — 
exp(-iHt) acts on coherent states as 

U(t)\z) = exp(-iai/2)| exp(-iai)z). (95) 

We can see that the correspondence between U(t) and V(t) is not 1- 
to-1, but 2-to-l. For example, suppose that t = 1 and a changes to 
a + 2kir. V(t) is unchanged by this transformation, but U(t) changes 



to (— l) k U(t). This is consistent with the relation (93). 

In general, the Maslov index changes if the gauge transforma- 
tion is topologically non-trivial. (This is considered to be a global 
anomaly. 0, |Tl|) The gauge transformation S(t) is an element of the 
fundamental group of Sp(2n, R), and this group is isomorphic to Z 
Hi, namely, 

ni(Sp(2n,R)) ~ Z. (96) 

Therefore each gauge transformation has an integer k as a winding 
number, and fi is changed by this transformation as /i — > fi + 2k. 
Hence we should use only topologically trivial gauge transformations 
when we transform A into a normal form. 

In the following subsections, we formulate canonical transforma- 
tions for quadratic path i ntegr als rigorously, and show normal forms of 



these path integrals. In § 3.3.1 , we introduce metaplectic group follow- 
ing the reference p2| . Metaplectic group is a unitary representation 
of symplectic group, and the 2-to-l correspondence between the two 
groups is the origin of the anomaly as we saw in the Harmonic oscil- 



lator model. In §3.3.2, we formulate canonical transformation using 



metaplectic gruop. In § |3.4| , we discuss the classical time evolution op- 
erator V(t). Since this is a symplectic matrix with a parameter t, we 
can regard V as a curve on the group manifold of Sp(2n, R). Canon- 
ical transformations are represented by deformations of this curve. In 



3. J: , we define normal forms of quadratic path integrals, and calculate 



them. Details of the calculations are shown in appendix M, and the 



results are summarized in §3.6. We derive the formula for the Maslov 



index of the n-repetition of a orbit in §3.7 



3.3 Canonical transformations for quadratic path 
integrals 

3.3.1 Metaplectic group 

In this section, we introduce Metaplectic group, which is a quantum 



counterpart of symplectic group 12 . Let us consider a linear canonical 



transformation specified by a symplectic matrix S £ Sp(2n, R): 

x = Sx. (97) 
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The metaplectic operator corresponding to S is the unitary operator 
M(S) specified by the relation 

M\S)xM{S) = Sx. (98) 

Metaplectic group Mp(2n, R) is composed of such unitary operators, 
and locally isomorphic to Sp(2n, R). Let us consider an infinitesimal 
canonical transformation 

S(e) = 1 - eA. (99) 

A satisfies 

A T J + JA = 0, (100) 

which can be written as 

A = Jh, (101) 
where h is a symmetric matrix. Corresponding metaplectic operator is 

M{S{e)) = 1 - l -^x T hx. (102) 



It is easy to verify that (Jl02J) and (|99j) satisfy (|98|), using the commu- 
tation relation f] 

[x a ,Xf}] = -iJ a p- (103) 
Let us consider a symplectic matrix V(t) specified by the equation 

jV{t) = -Jh(t)V(t), (104) 
V(0) = 1. (105) 

V(t) can be represented as 

V(t) = P cxp[- / A(t')dt], (106) 
Jo 

where 

A{t) = Jh{t). (107) 

The corresponding metaplectic operator U (t) is uniquely determined 
by the following equations: 

lj-U(t) = -iH(t)U(t), (108) 
at 



2 Note that we assume h — 1 in this section. If Ti ^ 1, [x a ,xp\ = —iTiJ a $ and M(e) 
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Hamiltonian 


eq. ol motion 


time evolution operator 


classical 


Hit) = \x T h(t)x 


J± x = h(t)x 
(Hamilton eq.) 


V(t) =Pexp[- /J A{i)dt] 
A(t) = Jh{t),V{t) e Sp{2n,R) 


quantum 


H{t) = ±x T h{t)x 


i&\<p) = H(t)\<p) 
(Shrodinger eq.) 


U{t) = Pexp[-i f*H(t')dt'] 
U(t) e Mp(2n 7 R) 



Table 2: Correspondence between classical and quantum systems for 
quadratic Hamiltonian. 



U(0) = 1, (109) 
H(t) = ^x T h(t)x. (110) 
U{t) can be represented as 

U(t) = P exp[-i / H(t)dt]. (Ill) 
Jo 

V(t) is the classical time evolution operator for the quadratic Hamil- 
tonian H(t) = x T h(t)x, and U(t) is the corresponding quantum time 
evolution operator. V(l) (classical time evolution operator for the 
period) is called the monodromy matrix. We summarize the corre- 
spondence between classical and quantum systems in Table |^. 

Although Mp(2n, R) and Sp(2n, R) are isomorphic locally, they 
are not isomorphic globally; there is a 2-to-l correspondence between 
them. Therefore if a closed curve on symplectic group 5(i)(5(0) = 
5(1)) has winding number k, the metaplectic operators corresponding 
to 5(0) and 5(1) are not necessarily the same, but 

M(5(l)) = (-l) fc M(5(0)). (112) 

This relation is the origin of the anomaly. Note that the relation be- 
tween Sp(2n, R) and Mp(2n, R) is very similar to the well-known rela- 
tion between 50(3) and SU(2). The point is that both Sp(2n, R) and 
50(3) are not simply connected and hence allow two- valued represen- 
tations. 



3.3.2 Canonical transformation 



Let us consider canonical transformations of quadratic path integrals 
using Metaplectic operators. In continuum notation, a partition func- 
tion 



Z = I Vx exp 



A(t) x(t)dt 



(113) 
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seems invariant under the canonical (gauge) transformation by {S(t)} € 
m(Sp(2n,R)): 

x(t) -> S(t)x(t), (114) 

A(t) - 5(t)A(t)5- 1 (t) + — (i), (115) 

5(0) = (116) 

However, this is not correct. We must go back to the definition of the 
path integral to make the discussion rigorous. 

Z is the trace of the time evolution operator U(t — 1), which is 
represented as a product of infinitesimal time evolution operators: 

^= i )=n( i -^(*i))+°(^) I ( ii7 ) 



where tj = 1/N. We obtain normal path integral by inserting the 
complete set 

1 = J d qj \ qj ){ qj \ (118) 

to ( |TT7| ): 

r N ' / i 

z = j W'Ui, <,, i -iiu,u h , • o (-^-) • rm» 



To formulate the canonical transformation corresponding to S(t), we 
insert 

1 = / dq 3 M(Sfe)) t |Q i )(9 i |M(5fe)) (120) 



instead of (118). Then (119) is modified as 



(121) 



where M 3 denotes M(S(tj)). 
S(tj-i) can be expanded as 



(ijJSfel-'UfeJ+O^V (122) 



Therefore 



M i Mj_ 1 = l-^|j5(t J -)^(* i )} ^+o(ij), (123) 
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M , //•:/, ;.W * , = MjAitJMj + ofjjY (124) 



= ^Kfid^MjX^MjipM} +o(Jj), (125) 
= Ih^it^S-^t^Sptitj^xs + O (Jj) (.126) 

Here we used M(S)xM(S)1 = S^x and ^S- 1 = S^-. Hence, for 

2<j<N, 

I 1 - jf 6 ^} M}_ x \q^) = (g.|l-A^ fe) | q ._ l)+0 (± 

(127) 

2<f(i) = x T k{t)x, (128) 

fc(t) = (S(i)-TMW)- 1 - J5(t)-^-(t). (129) 

This is the same as the classical canonical transformation. However, 
for j = 1, 



M\Mfo = (-l) fc + O I - 1 , (130) 

where k is the winding number of the closed path {S(t)}. Therefore Z 
is represented as 

N 

I I \ 

' . \ \ . \ ; 

JV 



z = /n^i<9ii i -^(*i)i9i-i)+°(^). ( m ) 

= (-!)"/ n d 9i<9il 1 -^fe)l«J-i>+°(^)'( 132 ) 

and the path integral (113) is transformed as a functional of A like 

ZL4 s (t)] = {-l) k Z[A{t% (133) 



where A 5 is the gauge field transformed by S. (133) means that the 
Maslov index is transformed as 

fi = n + 2k (mod4), (134) 

where (x is the Maslov index corresponding to A s (t). We show the 
stronger result 

\l = [i + 2k, (135) 

in the following subsection. 



3 Note that this simple correspondence between classical and quantum canonical trans- 
formation is due to Weyl ordering of the operators. If the Hamiltonian is written in other 
order, additional terms appear and the simple correspondence is lost. See, for example, 
chapter 3 of the reference 
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definition 


gauge transformation 


Hamiltonian 


H(p,q,t) = x T h(t)x 
h: symmetric matrix 


(S(t)-YHt)S(t)-i-JS(t)^(t) 


gauge field 


A(t) = Jh(t) 
A T J + JA = 


S(t)A(t)S(t)-i+S(t)^(t) 


curve on Sp(2n, R) 


(d/dt + A(t))V(t) = 0, V(0) = 1 
V(t) T JV{t) = J 


S(t)V(t)S(0)-i 



Table 3: Hamiltonian H, gauge field A and curve V. Those three have the 
same informations. S(t) belongs to Sp(2n,R) and satisfies S(l) = 5(0). V 
and S is continuous, but we don't require smoothness of V and S in this 
paper. Therefore H and A allow discontinuities. 



3.4 Curves on Sp(2n, R) 

To investigate quadratic path integrals in detail, let us examine the 
classical time evolution operator V, which is considered to be a curve 
on the group manifold of Sp(2n,R). {V(t)} is the curve which starts 
from the origin, and its end point is the monodromy matrix. V is 
defined by the equations 

pV(t) = j(±+A(t))v(t)=0, (136) 



dt 

V(0) = 1. (137) 
Note that the correspondence between V and If) is 1-to-l. V is deter- 



mined from If by ( J13q ) and (137), and Jp is determined from V by 
the relation 

A(t) = -^(t)V(t)-\ (138) 



See also Table |. 



3.4.1 Gauge transformation of V 

Let us consider how this curve is transformed by gauge transforma- 
tions. If the gauge field is transformed as 

A(t) - SWAQS®- 1 + 5(t)^-(t), (139) 



solutions of (136) are S(t)V(t) x c onstant. The constant should be 
chosen to be S^O)" 1 to satisfy (137). Thus V is transformed as 

V(t) -» SffiVtySiQ)- 1 . (140) 
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The gauge transformation group G is considered to be a set of 
closed curves on Sp(2n, R): 

G = {S\S : [0, 1] -> Sp(2n, R), 5(0) = S(l)}. (141) 

G can be decomposed into two parts, 

G = Lx N, (142) 

where L and N are subgroups of G defined as 

L = {l\l € G, 1(0) = 1}, (143) 

TV = {n\n : time independent element in G} (144) 

iV is isomorphic to Sp(2n, R). (142) means any S <G G can be written 
as 

S(t) = l(t)n, (145) 
where I € L and n £ N . I and n are explicitly written as 

n = S(0), (146) 

l(t) = SffiStO)- 1 . (147) 
The gauge transformation by I € L is written as 

V(t) - Z(t)V(t). (148) 

The monodromy matrix M(= V(l)) is unchanged by this transforma- 
tion. The transformation by n G -/V is 

V(t) — > nV(t)n _1 , (149) 

and M is transformed as nMn -1 . 

The elements in G and L are classified by winding numbers. We 
refer the subset of G (L) which has winding number k as G& (£*). Go 
(Xo) is a subgroup of G (L), and TV is a subgroup of Go- It is obvious 
that 

G = L x AT. (150) 

It is worth pointing out, in passing, that the time translation t — > 
t + <o is also regarded as a gauge transformation. This transformation 
changes V as 

™_J ^(i + io^o)" 1 (0<*<l-*o) ri - n 

1 n* + «o - 1) ^(to)- 1 a - to < t < i) ( j 

This is a gauge transformation given by 

om_J nt + toW)- 1 (0<t<l-to) Mt . 9 x 
6W ~ \ V(t + to - l)^(t)- 1 (1 - *o < * < 1) ( J 

which can be shrunk to a point continuously as to — * 0. Therefore the 
time translation is regarded as a gauge transformation in Go- 
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3.4.2 Continuous deformation of V and its Maslov in- 
dex 

Let us consider deformations of V. If the deformation is a gauge trans- 
formation, the monodromy matrix M is transformed as nMn -1 by a 
symplectic matrix n. The converse also holds. Let Vi and V2 be the 
curves which start from the origin. If there is a symplectic matrix n 
which satisfies V\(l) = riV^l)™ -1 , the matrix S(t) defined by 

S(t) = V 1 {t)nV 2 {ty 1 (153) 

satisfies 5(0) = 5(1) = n and Vi(t) = S(t)V 2 (t)S(0)- 1 . Therefore 
V\ and V2 are connected by the gauge transformation S. Thus we 
conclude that two curves are connected by a gauge transformation if 
and only if the end-points of two curves are in the same conjugacy class. 
As a special case, if two curves have the same end-point (monodromy 
matrix), there is a gauge transformation which connects two curves. 
In this case, n = 1 and S(t) = V\{t)V2{t)~ 1 is an element of L. 

Since there is a 1-to-l correspondence between lp and V, If) and 
its spectrum changes continuously as we deform V continuously. We 
can see relative changes of Maslov index by analyzing the points where 
the signs of diagonal elements of Jp changes. Q lp can be diagonalized 
explicitly, and what is important is that the diagonal elements are 
unchanged by gauge transformations. (See appendix |a|.) Let V\ and 
V2 be curves which are connected by a gauge transformation S. If 
S belongs to Go, S can be shrunk to 1 continuously. Therefore Vi 
and V2 can also be united continuously, and diagonal elements of lp 
is unchanged throughout this process. Thus we conclude that if two 
curves are connected by Go, two curves have the same Maslov index. 
However, If two curves are connected by Gk (fc ^ 0), two curves have 
different Maslov indices. Now we can prove the following important 
statement: 

If two curves V\ and V2 satisfy the relation 

V 1 (t) = S(t)V 2 (t)S(0)- 1 , (154) 

S G G k , (155) 
the corresponding Maslov indices \i\ and \x 2 satisfy 

Ai2=^i + 2fc. (156) 

4 Strictly speaking, we should determine the Maslov index from discretized matrix Jp N . 
However, when we deform V (or A) continuously, relative changes of the Maslov index 
are dominated by the low- frequency part of It) N , which is approximated well by lp UN is 
large enough. Therefore we can determine the relative changes of the Maslov index from 

JP- 
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Figure 3: Deformations of curves on Sp(2n,R). 



[proof] 

V2 can be deformed to V 2 , which is defined as 

Vz(t) = S(0)V2(t)S(0)- 1 . (157) 

Since the constant matrix 5(0) belongs to Go, Vi and V 2 have the 
same Maslov index. V\ and V 2 have the same end-point. (See Fig. ||) 



We already know that (156) holds for some cases. For example, let us 



consider the Hamiltonian 

H{p,q,a) = ~(p\+ql). (158) 

Let W\ and W 2 be the curves corresponding to H (p, q, a) and H (p, q, a+ 
2kn). Wi and W2 have the same end-point, and the difference between 



the Maslov indices of W\ and W% is 2k. (See section 3.2 and appendix 
|b|.) Then deform V\ continuously to V 2 through W\ and W2 (Fig. ||). 
The difference of the Maslov index between V\ and V 2 is the same as 
that of W\ and W2, because the change of the Maslov index during the 
deformation V\ — > W2 is canceled by that of W% — > V 2 if the end-point 
go through the same path during the two processes. Therefore the 
difference of the Maslov index between V\ and V2 is also 2k. (Q.E.D.) 



(156) means that two curves considered to be equivalent if two 
curves are connected by Go (not G) . We classify quadratic Hamiltoni- 
ans by this equivalence relation in the next subsection. 
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3.4.3 Universal covering of Sp(2n, R) and generalized 
monodromy matrices 

Let us consider the set of curves on Sp(2n, R) which starts from the 
origin. We refer this set as F: 

F = {V\V : [0, 1] -> Sp(2n, R), V{0) = 1}, (159) 

which is isomorphic to the set of quadratic Hamiltonians and the set of 
gauge fields. Two elements V\ and V 2 in F is considered to be equiva- 
lent if they are connected by a gauge transformation which belongs to 
G . We express this relation as 

Vi ~ V 2 . (160) 

G 

If Vi ~ V 2 , Vi and V 2 have the same Maslov index. Therefore the 
Maslov index is determined uniquely on the quotient set F / ~° . 

Before turning to the division by Go, let us discuss the division 
by L , which is a subgroup of G . If V\ and V 2 are connected by 

an element in Lq (i.e. V± ~ V 2 ), two curves have the same terminal 
points and the closed curve made by these two curves can be shrunk 
to a point continuously. The converse also holds. (It is obvious that 

{Vi(t)V 2 (t)~ 1 } belongs to L .) Therefore the quotient set Fj ~ is 
the set of homotopy classes, which is simply connected by definition. 
Therefore this is the universal covering space of Sp(2n, R): 

= Sp(2n,R). (161) 

An element in Sp(2n, R) is specified by a symplectic matrix (mon- 
odromy matrix M) and an integer (winding number k). (See Fig. ||) 
This pair M — (M, k) is regarded as a generalized monodromy matrix. 
We give a definition of the winding number in §3.5.2|. 



Let us discuss the group structure of F and Sp(2n,R). We can 
define the product of two elements V and W in F as 

{VW)(t) = V{t)W(t) (0 < t < 1). (162) 

It is obvious that F becomes a group by this product. Therefore we 
define the product in Sp{2n, R) = Fj ~ as 

[V][W] = [VW], (163) 

where [V],[W] and [VW] denote the homotopy classes of V,W and VW. 
It is easy to verify that this product is well-defined, and Sp(2n, R) also 
becomes a group by this product. 
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We can define another product in F as 
(V* W)(t) - 



W(2t) (0 < t < 1/2), 

V(2t-l)W(l) (1/2 < t < 1). 



(164) 



Note that this product doesn't satisfy the associative law, i.e. V * 
(W * X) 7^ (V * W) * X. However, VW and V * W is homotopic. 
We can define a map / : [0, 1] x [0, 1] — > Sp(2n, R) which satisfies 
f(t,0) = (VW)(t) and f(t, 1) = (V * W)(t) as 



f{t, s) = V(t, s)W(t, s) (0 < t < 1, < s < 1), 
f 1 (0<t<s/2), 



V ^Ef (*/2<*<l). 



W(M) = ( (0<*<l-/2), 
1 W(l) (l- S /2<t< 1), 

Therefore, in Sp(2n,R), 

[V][W] = [VW] = [V*W}. 

3.5 Normal forms of quadratic Hamiltonians 

It was observed in the preceding subsection that 

L 



(165) 
(166) 

(167) 
(168) 



F/ 
F/ 



Ln 



= Sp(2n,R), 
= Sp(2n,R). 



(169) 
(170) 



G Go 

We are now ready to consider Fj ~ and Fj ~ . Since G 



Lx N, 

Fj H is the set of conjugacy classes of Sp(2n, R). Therefore an element 

in Fj is specified by a conjugacy class of Sp(2n, R) and a winding 
number. In the following, we show representatives of these quotient 
sets. 

In this subsection, we change the definition of the 2n x 2n matrix 
J as 

( j \ 

j 



J = 



3 / 



3 = 



1 

-1 



(171) 



(172) 



for convenience. 
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3.5.1 Normal forms of monodromy matrices 

First we consider representatives of the quotient set 

F/ ~ = {conjugacy class of Sp(2n, R)}. 



(173) 



Our task here is to transform a given monodromy matrix M G Sp(2n, R) 
into a simple normal form SMS^ 1 , where S is also a symplectic ma- 
trix. In other words, we choose the special symplectic basis in which 
M has a simple form. 

Let us consider the eigenvalue equation 



Mx = Xx. 



(174) 



First we assume that M has no degenerate eigenvalue. If A is an 
eigenvalue of M, A* and 1/A are also eigenvalues of M [l3j] . (Fig.||) 

Therefore we can classify eigenvalues of M into four types, and M 
become a block diagonal matrix in the basis made of the eigen vectors. 
(If some of the eigenvectors are complex, we use the real parts and 
imaginary parts of them as the basis.) 



M = 



V 



mi 



m 2 



\ 



mi J 



(175) 



In the following, we list up the four types of eigenvalues and the 
normal forms of the monodromy matrix. 
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1. elliptic type: A 



= e ±la . 



(0<a < 2tt) 



cos a 
sin a 



2. hyperbolic type: A = e ±/3 . 



m 



— sm a 
cos a 



3. inverse hyperbolic type: A 



,±0 



m 



(176) 



(177) 



(178) 



4. loxodromic type: A = e ±M±/3 . (0 < a < n) 



V 



e 13 cos a 



sin a 



-e' 3 sin a 



e ^ cos a 


-e - ^ sin a 



e ^ sin a 



cos a 



(179) 



Here, a and /3 are real numbers. 

If M have degenerate eigenvalues, we must treat such cases sepa- 
rately. In this paper, we treat only the most important case, parabolic 
type. In this case, the eigenvalue 1 is doubly degenerate and we can 
choose the basis x a (a = 1, 2) as Fl 



Mx\ = xi, 
Mx 2 = x 2 - 

The monodromy matrix of this part is 



JXi. 



1 -7 
1 



(180) 
(181) 



(182) 



In the elliptic case, the normal forms corresponding to a and 2tt — a are not conjugate, 
though two matrices have the same eigenvalues. However, in the loxodromic case, two 
matrices corresponding to a and 2tt — a are conjugate. That's why we choose a as 
< a < ty in this case. See Appendix A. 2.1 and the reference [la], §42. 

6 7 can be transformed to ±1. {x\ —* v/py|a;i, X2 — » ( 1 / \/rTl ) 33 2 . ) However, we leave 
7 as a parameter for convenience. 
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3.5.2 Normal forms of Hamiltonians 

Now we can choose representatives of Fj ~ and Fj ~ . Normal forms 
of Hamiltonians corresponding to them arc also determined. 

First we show the representatives of Fj ~ and the normal form of 
the Hamiltonian corresponding to the types of of eigenvalues. 

• elliptic type 

cos at — sin at 
sin at cos at 



V(t) 



(183) 
(184) 



• hyperbolic type 



V(t) 



oft 

e" 
H = -f3pq. 







at 



• inverse hyperbolic type 
V(t) = 



cos 2Tit — sin 2nt 

sin 27rf cos 27rf 

e /3(2t-l) 

_ e -0(2t-i) 



(0 < t < 1/2), 
(1/2 < t < 1) 



| 71" (P ; 



2 + 9 2 ) (0<t<l/2), 
2/3pq (1/2 < i < 1). 



(185) 
(186) 

(187) 
(188) 



• loxodromic type 

/ e^* cos at 

V(t) = 



V 



e' 3 ' sin erf \ 

e _/3t cosat e - ' 3 ' sin erf 

—e^* sin at e^* cos at 

-e _/3t sinat e~^*cosai J 

(189) 

= a(pi<7 2 - P2?i) - + P292)- (190) 



• parabolic type 



K(t) 



1 



2 y 



-7* 
1 



(191) 
(192) 
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Figure 6: Normal form of V 



Let us consider a given time-dependent quadratic Hamiltonian H 
and the corresponding curve V. We can transform V into Vb, which is 
a direct sum of the normal forms, using the gauge transformation in 
G: 

V (t) = S(i)V(t)S(0)- 1 . (193) 
Then we can choose a normal form of V as 



Vk = Vo* U k . 



(194) 



Here, k denotes the winding number of S, and Uk is a closed curve 
which start from the origin and have the winding number k. (See 
Fig.0.) We can choose a normal form of Uk as, for example, 



U k (t) 



t 



\ 



cos 2nkt 
sin 2iikt 



— sin 2nkt 
cos 2nkt 



\ 



1 / 



We use this normal form Vq * Uk as a representative of F/ 
corresponding Hamiltonian is 

H(r,nt\-l 2 Mp? + 9?) (0<t<l/2), 

lP ' 9 ' ' ~ 1 ffo(p, 9, 2* - 1) (1/2 < t < 1) 
Here, £/o is the Hamiltonian corresponding to Vq. 



(195) 



The 



(196) 
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stability 


A 


diagonal elements 


|Dct^> | 


index 


elliptic 


e ±ia 


—a + 2rn7r 


4 sin 2 a/2 


1 


hyperbolic 


e ±P 


±f3,±^/(3 2 + (2rmr) 2 


4sinh 2 /3/2 





inverse hyperbolic 


_ e ±/3 


± v / f3 2 + (2m + 1) 2 tt 2 


4 cosh 2 p/2 


1 


loxodromic 


e ±ia±/3 


± v / f3 2 + (a + 2rmr) 2 


4(cosh/? — cos a) 2 





parabolic 


1 


0,-7, -7/2 ± V(7/2) 2 + (2m7r)2 


(Dct p = -7) 


(l/2)sgn 7 



Table 4: The absolute values of the functional determinants, diagonal ele- 
ments of If) , and the Maslov indices for the normal forms. m runs over all 
integers in elliptic and loxodromic cases. In hyperbolic and parabolic cases, 
m > 1, and in inverse hyperbolic case, m > 0. 



3.6 Absolute values and Maslov indices of quadratic 
path integrals 

We can reduce a curve V to the normal form Vq * Uk by means of a 
gauge transformation in Go . The absolute value and the phase (Maslov 
index) of the path integral are the invariants of this transformation, 
and we can easily calculate the path integral of this normal form. 

The absolute value of the path integral corresponding to Vo*Uk is 
the same as that of Vq, and the Maslov index of Vo*Uk is (index of Vq) + 



2k. Since Vq is a direct sum of the normal forms shown in 3.5.2, the 
calculation of the path integral is reduced to the calculation of the 
path integrals of the normal forms. We summarize the results of the 
calculations in Table |^. For details of the calculations, see appendix 
§ 

The functional determinant corresponding to Vq is the product of 
the functional determinant corresponding to the normal forms, and 
the Maslov index of Vq is the sum of the Maslov index of the normal 
forms and 2k. Therefore, if the monodromy matrix has no degenerate 
eigenvalue, we obtain the formula 

|Det#>| = |det(M-/)|, (197) 

H=p + q + 2k, (198) 

where p (q) is the number of elliptic (inverse hyperbolic) pairs of eigen- 
values. The partition function become 

Z = J Vxexp ^x T Ipx , (199) 

e -*f P 

(200) 



V|det(M- J)| 
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If the monodromy matrix has parabolic blocks, we must add J2i 5 s g n 7i 
to the Maslov index. In this case, Ip has zero-modes, and Det^Z> = 
0. In the space where the zero-modes are removed, the functional 
determinant become 



Det^|=V|det(M'-J)l | 7i |. (201) 




Here, Det denotes the functional determinant in the space without 
zero-modes, and M is the monodromy matrix whose parabolic parts 
are removed. Integration with respect to zero-modes will be discussed 
in §0. 



3.7 Repetitions of an orbit 

Let us discuss repetitions of an orbit. This is the problem of calculating 
the partition function 

Z(n) = Tr U(T = n) = Tr{U(T = 1)™}. (202) 

If the classical time evolution operator corresponding to Z(l) is V, 
that corresponding to Z(n) is V n — V * V *....* V. Let the normal 
form of V be Vq * Uk- Then we can transform (Vq * Uk) n into V n * Ukn 
continuously. (See Fig. |?].) Therefore the problem is reduced to the 
calculation of the repetition of the normal forms. 

The Maslov index of the normal forms of hyperbolic, loxodromic 
and parabolic type are unchanged by the repetition. For elliptic type, 
the winding number increases when na goes over 2tt, where n is the 
number of the repetition. Therefore the Maslov index of this part is 
1 + [na /2tt] . The normal form of inverse hyperbolic type can be written 
as Vh * Ui , where Vh is the normal form of hyperbolic type. Therefore 
n-repetition of this normal form can be transformed to V£ * Ua . If n 
is even, this is hyperbolic, and if n is odd, this is inverse hyperbolic 
type. The Maslov index of this part is n. 

Thus we obtain the formula for the Maslov index of the n-repetition: 



p / 

i=l 



2vr 



qn + 2nk. (203) 



Here, at denotes the stability angle of i-th elliptic block. Note that 
jj, n ^ n/ii if the orbit has elliptic blocks. 
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Figure 7: (Vq * Uk) n can be transformed into Vq * U-zkn- 
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4 Trace of the resolvent operator 



In this section, we discuss the trace of the resolvent operator 

9(E) = Tr Fbr 



(204) 



to discuss eigenvalues of the time-independent Hamiltonian H. This is 
obtained by Fourier transformation of the partition function: 



9(E) 



1 



e iET ' n Z{T). 



(205) 



We can find eigenvalues of H from the poles of g(E). 

Z{T) become a sum over periodic orbits, as discussed in the pre- 
vious sections. In this case, each orbit has at least one zero-mode 
corresponding to time-translation symmetry. Therefore we first dis- 



cuss the integration with respec t to zero- modes in §4T. Then we carry 
out Fourier transformation in §4^2, a nd fi nally we give a formula for 
the sum over repetitions of a orbit in §4.3. 



4.1 Integration with respect to zero- modes 

If the system has continuous symmetries, jp has zero-modes and we 
must execute the integration with respect to them. In our formalism 
based on the path integral representation of the partition function, it is 
easy to treat this integration because the integrand (— i det If) ) -1 / 2 exp [iR\ 
is a constant. Therefore the contribution from a continuous family of 
orbits is 

V 

^=^= exp 



-zdet p 



i 



(206) 



where V is a volume factor of the family of orbits. 

Hereafter, we consider an isolated orbits in a time-independent 
Hamiltonian system. In this case, the orbit has a parabolic block 
corresponding to the time translation symmetry. 

The parabolic block become 



(207) 



1 

dT -, 
" dE 1 



7 Note that the notation of the vector in phase space is different from M. In this paper, 
we note p (or E in this case) first, and q (or t) second. This is opposite of the notation in 
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which is conjugate to 



1 d/r 

dE 



1 

Therefore the contribution from this block is 

T 



(208) 



(209) 



(See appendix ||.) and the amplitude factor of this orbit is 

K = 



2nh^y/\det{M' -I)\ 



(210) 



Here, M denotes the monodromy matrix whose parabolic part is re- 
moved, and the Maslov index fi become 



1 



H=p + q + 2k- -sgn — 



dT 



dE 



(211) 



4.2 Conversion from T to E 



Let us discuss the conversion from a fixed period T to a fixed energy 
E. The partition function for a fixed T is approximated by a sum over 
periodic orbits as 



Z(T) = K exp ~ j> pdq - Hdt 



(212) 



where the amplitude factor K is given in (21C). The trace of the 
resolvent g(E) is the Fourier transformation of Z(T): 



(213) 



g(E) = / dT if exp 



-{ET+cp {pdq - Hdt) 



Let us evaluate this integral by the stationary phase approximation. 
The stationary phase condition 



d 

~dT 



ET + j>(pdq - Hdt)\ = 



leads to 



H(p(t),q(t))=E. 
The second derivative of the phase is 



— \ET+ j> (pdq - Hdt) 



1 dE 



(214) 
(215) 

(216) 
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and the contribution from this factor cancel the contribution from the 
parabolic block (207). Therefore g is approximated as 



g(E)=J2 K ' exp 



S = i> pdq, 



K 



(217) 

(218) 

(219) 
(220) 



V|det(M'-7)|' 
\i = p + q + 2k. 
This is the well-known Gutzwiller's trace formula. 

4.3 Sum over repetitions 
4.3.1 One dimensional case 

Let us discuss the sum over repetitions of a orbit using the formula 
(203). First, we consider the one dimensional case. In this case, there 
is only one type of periodic orbit. The monodromy Matrix has only 
one block corresponds to time translation symmetry, and the winding 
number of the orbit is 1. (See Fig. ||.) Therefore the Maslov index of 
the prime periodic orbit is 2 and the index of the n-repetition is 2n. 
The trace of the resolvent g is approximated as 



9(E) 



9o(E) 
9o(E) 



OO 

^Texp 

n=l 



-S(E) 



T 



expi[S{E)/h-n] 
l-expi[S(E)/h-%]' 



(221) 
(222) 



Here, go(E) is the Thomas- Fermi term which represents the contribu- 
tion of the zero-length orbit. T and S are the period and the action of 
the prime periodic orbit. 

Since the poles of g[E) represent the eigenvalues of the Hamilto- 
nian, we obtain the Bohr-Sommerfeld quantization condition: 



S{E) = <j)pdq= (2n + l)7rft. 



(223) 
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Figure 8: The periodic orbit of the one-dimensional system. The local 
coordinate around the orbit rotate once in the phase space when the orbit 
goes around the course once. 



4.3.2 Elliptic orbits 

Let us consider the periodic orbits where all eigenvalues of the mon- 
odromy matrix is elliptic except the parabolic block corresponds to the 
time-translation symmetry. In this case, 



|det(M' --01 = II | 2sin ' 



(224) 



where p is the number of elliptic blocks, and ay is the stability angle 
of the j-th elliptic block. The Maslov index is 



fi=P + 2k, 



(225) 



where k is the winding number. Hence the amplitude factor of this 
orbit is 

(— l) fc T 



K rr^sinf 

From the formula fl203| ), we obtain 

p 

m» = Yl ( 



i 



na, 



2tt 



2kn 



and 



K„ 



{-l) kn T 
IlLi2isiii 



na i 
2 



(226) 



(227) 



(228) 
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for n- repetition of this orbit. 

Therefore the sum over repetitions of this orbit is 



^2 K n exp 



n=l 



-nS{E) 



p oo 



(mj + ±)nctj 



exp 



n-1 



j — 1 mj—0 



^nS(E) [229) 



n=l { mj } 



< /. - • ^ | ;;/ - - ; n , 



^E 



exp • 



S(£?)-{Awr + E?=i to + ?)<*,}*] 



/ F 7 =j 

{mj} 1 - exp * [s(25) - {/c^ + £* =1 (mj + i) ajj h] 



Here, we used 



sin § ^ 



(232) 



2 m=0 

Therefore the quantization condition is 



S(E) = | (2m + k)n + ^ (m^ + 1^ ttj J ft, (233) 

where n is an arbitrary integer and m,j is an arbitrary integer which is 
not negative. Note that the quantization condition depends on if the 
winding number k is even or odd. 



4.3.3 Unstable orbits 

Let the monodromy matrix of a periodic orbit have p elliptic blocks, 
q hyperbolic blocks, r inverse hyperbolic blocks, s loxodromic blocks 
and a parabolic block corresponding to the time translation symme- 
try. Then the eigenvalues of the monodromy matrix can be written as 
exp(±z£) using the stability angle £, and £ is represented as 



Ceje = a eJe {l<j e <p), 

Ch,j h = i/3h,j h (i < jh < q), 

(ih,j ih = iPih,j ih + 7T (1 < j ih < r), 

C,l,h,± = ±a i,n + , (! < 3l < s ) 

Cp = 1. 



(234) 
(235) 
(236) 
(237) 
(238) 



Here, < ct e ,j e < 2tt and < aij l < it. (3h,j h , Pih,j ih and are 
positive. p + q + r + 2s + l = n where 2n is the dimension of the phase 
space. 
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The amplitude factors of this orbit and its repetitions are 



Kn 



(-1) 



knr 



l,n i 



where 



K, 



1 



n 



2i sin : 
1 



K, 



h.n 



n « =i 2sinh2^ 

IX=i 2icosh ~^ 

(-1)" /2 



n; =1 



2 sinh ■ 



(n: odd) 
(n: even) 



K, 



1 



(239) 

(240) 
(241) 

(242) 
(243) 



nj=i 2(coshn/3; J - i — cos nay,) 

The sum over repetitions of this orbit can be calculated in the same 
way as elliptic orbits, and the result is 



^2 K n exp 



n=l 



-nS(E) 



T 



m 1 



exp[j-S c (E, m)] 
exp[iS c (E,m)] 



where 



S c (E,m) = S(E) — < kir 



n-l 

E 

3=1 



3 ) o 



(244) 



(245) 



Here, m = (mi, rri2, m„_i) is the n— 1(= p + g + r + 2s) dimensional 
vector whose components are non-negative integers. Q represents one 
of stability angles in Eqs. ( ^34[ )-( p37| ). Therefore the quantization 
condition is 



S(E) = <^ (2m + fc)7r 



n-l 

E 



(246) 



where mo is an arbitrary integer. If some of stability angles have 
imaginary parts, RHS of ( 246 ) have imaginary parts and eigenvalues 
become complex numbers. 

In chaotic systems, most of the periodic orbits are unstable and 
the eigenvalues determined by the condition (246) also have imaginary 
parts. Such phenomena occur even if the Hamiltonian of the system is 
a Hermite operator and leads to false results. 

This means that, for unstable orbits, it is not enough to consider 
only one orbit and its repetitions. We should take into account the 
correlations of unstable orbits. This is a difficult but interesting future 
problem. 
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5 Conclusion 



We have investigated the semiclassical trace formula using phase space 
path integral. We derived the semiclassical trace formula by applying 
the stationary approximation to the phase space path integral of the 
partition function. 

We classified the quadratic path integrals around the periodic or- 
bit. This is equivalent to the classification of connections on the vector 
bundle over S 1 with the structure group Sp(2n, R). However, if the 
canonical transformation is topologically non-trivial, the Maslov index 
changes. Therefore we should regard two connections as equivalent 
only if they are connected by the topologically trivial canonical trans- 
formation. The results of the classification by this equivalence relation 
is shown in §3.5.2. 

We defined normal forms of time-dependent quadratic Hamiltoni- 
ans, and calculated quadratic path integrals. We also derived a formula 
for the Maslov index of the n-repetition of the orbit in §3.7, and derived 
the Bohr-Sommerfeld type quantization condition in §4.3. 

Throughout this paper, we concentrated on the quadratic path in- 
tegrals around periodic orbits. However, the analysis of higher order 
terms is often needed when we investigate realistic systems like nuclei 
and atomic clusters. Phase space path integrals used in this paper will 
supply powerful tools for such analysis. 
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A Calculations of diagonal elements in the 
continuum limit 



A.l Diagonalization of D 

First we consider the following eigenvalue equation: 



Dx = ( — + A(t) ) x = ex. 



(247) 



This equation is easy to solve because the time derivative has diagonal 
form. Let £^ be the eigenvector of the monodromy matrix M, and £ 
be the stability angle: 

M£ c = e-«€ C - (248) 
Then the solutions of (248) arc written as follows: 



(t) = exp(e c>m t)y(i)£ c , (249) 

Dx^.m = ec,mXQ^ m , (250) 
H,m = (C + 2m7r)i. (251) 
Here, V(t) is the classical time evolution operator which satisfies 



DV(t) = 0. 

We can calculate |Det#>| = |DetL>| from ( |25T| ): 



iDctDl 



n oo 



\{ \{ (2m7r + 0)(2m7r-0) 

l — l m— — oo 

71 OO 

n n ( 2 ^) 2 -c 

/=1 m— — oo 

no 2 n(2m 7 r) 2 (i 



1=1 m=l 



1=1 



n( 2 »4 



\2to7t 
= |det(M-7)| 



(252) 

(253) 
(254) 

(255) 
(256) 



Here, we used the formula (407) and ignored the diverging constant 
rim=i 2m7r. We can see from the calculations in Appendix [b| that this 
constant corresponds to 



N-l 



1 n „ . ^ m 

— TT 2 sin = 1 

N J- 1 , N 



(257) 



m— 1 



in the discrete formalism. 
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A. 2 Diagonalization of 



In this subsection, we diagonalize p using the basis given in A.l 



A.2.1 Elliptic type 

Eigenvalues of the monodromy matrix are 

Mx± = e Tia x ± (0 < a < 2tt), (258) 

(x+ = x*_). a is taken < a < tt if [x, Mx] > 0, and it < a < 2tt if 
[x, Mx] < 0, ^| where x is an arbitrary real vector in the space spanned 
by x + and x . See also [jL3| (§42). We rewrite x + as x + — x Re + ixi m 
where x Re and real vectors, then we obtain 

MxR e = x Re cos a + xi m sin a, (259) 
AIxj m = x Re sin a — xj m cos a. (260) 

The inequality 

[x Re ,x Im }>0 (261) 
follows from the definition of a. We use the normalization condition 

[x Re ,x Im ] = 1, (262) 

[x+,x-]=-2i. (263) 
The eigenfunctions of D are 

x±,m(t) = exp(e± im t)V(t)x±, (264) 

e ±,m = (±ct + 2rair)i. (265) 
We define the real functions x^m^Jm,™ as 

x Re , m = ~(x +tm + X-^ m ), (266) 

xi m . m = —(x + m -x-- m ) (n = 0,±1, ±2, ....). (267) 
2i 

These functions satisfy the normalization and the orthogonalization 
condition: 

/ dt [x Re (268) 
Jo 

/ dt [XRe, m ,X Re , n ] — I dt [xi m ,m, Xim,n] = 0. (269) 
JO JO 



3 The product [x,y] means x T Jy. 
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Matrix elements of Jp are 

/ dtxT Re. m J P X n 

Jo 



= -(a + 2mn)S m . n , (270) 

/ dtx^ mm ]p ccj m ,„ = -(a + 2mir)5 m ,n, (271) 
Jo 

I dtx^ e m Tp xi m ^ n = 0. 
Jo 



(272) 



Therefore Tp is already diagonalized by this basis. 
The functional determinant of this part is 



21 



2,73) 



i— — oo ^ m= 1 J m=l 

- 4 sin 2 -. (274) 



Here we used the formula (407). 
A. 2. 2 Hyperbolic type 

In this case, the eigenvectors of the monodromy matrix satisfy the 
equation 

Mx± = e ±/3 x ± (275) 
where (3 is a real number. The normalization condition is taken as 

[x+,x_] = l. (276) 

We define c± jTn ,x± tm (t) as before: 

a; ±!m = exp(e± l77l f)V(t)a;±, (277) 

e±,m = T/3 + 2mm. (278) 
We define real functions a; c .±,m, ^s,±,m as follows: 

■^c,±,m = ~ ^(«E±,m "1" «E±,— m)j (279) 
Sj,±,m = -r-7=(x±,m - X±- m ) (m = 1,2,...). (280) 

zy2 

Since £c± t o are already real, we don't have to redefine them. Matrix 
elements which have non-zero values are 



/ dtxl + m p: 
Jo 



= /3, (281) 
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JO 

I dt x s _j_ m lfl x c ^ rn 
Jo 

/ dtx +0 ipx-.o 
Jo 



—2mir, 

2TO7T, 



(282) 
(283) 
(284) 
(285) 



Therefore Tp is represented by the matrix d m in the space spanned by 



In the case m = 0, 









/* 


2rwr 








— 2m7r 


/? 


/5 


— 2m7r 








2m7r 


/? 









(to > 1). 



d 



/3 
/? 



The solutions of the eigenvalue equation det(d m — XI) = are 



A = ±v//3 2 + (2tott) 2 (m > 1), 
A = ±/3 (m = 0). 

The solutions for to > 1 are doubly degenerate. 



The functional determinant of this part is (using (408)) 

oo oo 

ndetd m = (3 2 J] {/3 2 + (2to7t) 2 } 2 , 

m— 1 

~ 4sinh 2/? 



m=0 



(286) 



(287) 



(288) 
(289) 



(290) 
(291) 



A. 2. 3 Inverse hyperbolic type 



In this case, eigenvalues of the monodromy matrix are negative real 
numbers: 

Mx± = -e ±p x±. (292) 
The normalization condition is taken as 



[x + , X—] = 1. 

Eigenvalues of D are 

e±,™ = T/3+ (2to+ 



(293) 
(294) 
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We define real vectors x c ± m , x s ± m as follows: 



f 



(295) 



a^s.i.m = T _ 7=( a; ±,m - aJ±,-m-i) (m = 0, 1, 2, ....) (296) 
zV2 

The matrix elements of If) can be calculated in the same way as the 
hyperbolic case. In the space spanned by x Cj ± j7n and x St ± >m , Jp is 
represented as 

/ (2m + 1)tt \ 

-(2m + l)n (3 

[3 -(2m+l)7r 

V (2to+1)tt (3 



Eigenvalues of c? m are 



A = ± v //3 2 + (2m + l) 2 7r 2 



(297) 



(298) 



Each solution is doubly degenerate. The functional determinant of this 
part is 



Y[d m = J] {/9 2 + (2m+l) 2 7r 2 } : 



m— 



m=0 



4 cosh 



(299) 
(300) 



Here we used (409). 

A. 2. 4 Loxodromic type 

Eigenvalues of the monodromy matrix are 

Mx ± . ± = e* ia± Px±,± (x±, ± = x* i± ). 
Normalization is taken as 

-2i 
-2i 

The eigenvectors of D are 

£C±,±,m(*) = exp(e± j ± >m i)V(*)a'±,±> 

where 

e±,±,m = T P + 2m7ri. 



(301) 

(302) 
(303) 

(304) 
(305) 
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We define real basis £c_R e ,±,m, xim±m as 



1 



These vectors are normalized as 

[^i?e,+,m) *^/m.-,n] Re, — ,m; '^'/m, + ,n] ^m. 

Other combinations are equal to zero. 

The matrix elements which have non-zero values are 



(306) 



dt x 



[ dtxj 
Jo 



i(x+,±, m - x_, ±i _ m ) (m = 0, ±1, ±2, ...).(307) 



(308) 

(309) 
(310) 
(311) 
(312) 



/ dt x\ e Jr m lpxi m ^-^ 
Jo 



-{a + 2to7t), 



= /3, 



-A 



m, + ,m 



} xim,-,m = -(a + 2m7r). 



Therefore ^ is represented in the space spanned by XR e .±.m, xim.±.> 
as 



/ -(a + 2TO7r) 

0-/3 
-(a + 2m7r) -/3 

\ (3 -(a + 2mir) 



/3 

-(a + 2m7r) 





Eigenvalues of d m are 



(313) 



(314) 



Since the determinant of d m is {/3 2 + (q + 2m7r) 2 } , the functional 
determinant of this part is 



rf 2 



Yl {(3 2 + {a + 2mn) 2 } 2 = \[ (a + 2mi) 4 f[ |l 

i— — oo m— — oo m— — oo ^ 

~ 4(cosh/3 — cos a) 2 . 



(a + 2m7r 



(316) 



Here, we used the formula (J407]) and (410) 
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A. 2. 5 Parabolic type 

If there are eigenvalues of the monodromy matrix which are equal to 
1, the operator Tj) has zero- modes. 

Here, we treat the most important case, the parabolic type. In this 
case, there are bases x\,x 2 which satisfy 

Mx\ = x 1: (317) 
Mx 2 = X2--JX!. (318) 

Here, the normalization is taken as 

[x u x 2 ] = l. (319) 

The monodromy matrix is represented in this space as 



1 - 7 
1 



We define a matrix B which satisfies 



(320) 



Bxi = 0, (321) 
Bx 2 = -jxi. (322) 

This is a generator of M . (e B = M .) 

We can diagonalize Ip using these bases. £Ci, m , x 2 , m arc defined as 

x a . m {t) = cxp{2mirit)V(t)aq>(-Bt)x a (a = 1,2). (323) 

D acts on these bases as 

Dxi_ m = 2mirixi. m , (324) 
Dx 2t , n = 2mirix 2 . m + jx ltm . (325) 

We define real bases x CAm ,a; Si ^ m as 

Xc,a.m ^^{x a .m ~t~ X a — m ) (326) 

Xs,a.m ~ 7={x a .m X a — m ) (?71 1, 2, ...). (327) 

x\fi and #2,0 are already real. 

Then the matrix element which is not equal zero are 



L 



i 

dtx c,l,mP X s,2,m = 2m7T, (328) 
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/ dtxJ lrn lpX c ,2. m = -277J7T, 

Jo 

/ dt X T %m Jp X Ci2 ,m = -1, 

Jo 

/ dtxj 2 Jpx s , 2 ,m = -7- 
Jo 



Therefore, for m > 1, 


















— 2m7r 





— 2m7r 


-7 


2rmr 











-7 / 



dct d m = (2to7t) 4 . 



Eigenvalues are 



+ (2tott) 2 . 



Each solution is doubly degenerate. 
For m = 0, 

d ° = (o -7 
The eigenvalues of this part is 



Since If) has a zero-mode, 



A = - 7 ,0. 



Dct^ = 0. 



(329) 
(330) 
(331) 



(332) 

(333) 
(334) 

(335) 
(336) 
(337) 



The functional determinant in the space where the zero-mode is re- 
moved is 

det'D = -7. (338) 
The zero-mode should be integrated separately. 
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B Direct calculations of some path inte- 
grals 

B.l Elliptic type 

Let us consider the normal form of the Hamiltonian of elliptic type: 

H=^(p 2 +q"). (339) 

The path integral for this Hamiltonian can be calculated directly. The 
eigenvalue equation of lp N is 

- jjPj + Qj -Qj-i = tPj, (340) 

Pj - Pj+i ~ (ij-i + 2 1j + Qj+i) = e 1j ■ ( 341 ) 

Since the Hamiltonian is time- independent, we assume that the solu- 
tions can be written as 

(S) = (") eXP (^) ^ = 0,l,..,iV-l). 042) 
Then the eigenvalue equation becomes 



(343) 



e satisfies the following quadratic equation: 

— 1 + cos 2 — e + — cos 2 4 sm 2 — = 0. (344) 

N V AT/ VAT/ N N V 7 



e 2 



We can find all eigenvalues by solving this equation. 
Let us calculate det Ip N . 



N-l 

detp N = 



u{^y^i-^i}- («») 

r=0 k ' 



^n(-^)n(i-^>«> 

r=l r— 1 v 



The first part of the product can be calculated using the formula (414): 

n>S-j£i- (347 » 
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For the second part, if N is large enough, 



N-l 



n u 



a 2 nr 



n 



r=l 



2 \ 2 

4N^ COt N 



r=l 



n 



We can calculate this infinite product using the formula (407): 



n i 



Therefore 



Sin 7TX 



lim | det Tp N | = 4 sin 2 — . 

AT— >oo 2 



,(348) 
(349) 

(350) 
(351) 



We can calculate the Maslov index using the quadratic equation 

2 „ irr 



(344). If 



iV 



4sm 2 — <0, 



(352) 



the solutions of the quadratic equation have opposite sign. Therefore 
this part doesn't contribute to the Maslov index. The problem is the 
case where 



/ a \ z o nr 9 nr 

— cos 2 4 sin — 

\nJ n n 



> 0. 



(353) 



In this case, the solutions of the quadratic equation have the same sign. 
If a is positive (negative), the coefficient of the linear term in ( |344| ) 
is also positive (negative) and both solutions are negative (positive). 
Therefore such pair contribute +1 (—1) to the Maslov index. 



(354) can be rewritten as 



N 



Tr) I" 4 - tan 2 - >0 (0 < r < N 1) 



N 



(354) 



r = always satis fies this inequality, and if r satis fies (354 ), N — r 
also satisfies ( (354|) . If N is large enough and r < N/2, (354) is reduced 
to 



1 



2nr 

a 



> 0. 



(355) 



The number of 0) satisfying this inequality is [^-] . Therefore, the 
Maslov index is 

li = X + 2 \^- \ . (356) 
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B.2 Hyperbolic type 

The normal form of the Hamiltonian of the hyperbolic type is 

H = -0pq. (357) 

This case can be treated in the similar way to the elliptic type. The 
eigenvalue equation is 

i-e-ity + JL(i + e -i*fr') } f p 



1 -e ia *r + JL(i + e 4 Tr) J\q 



(358) 



Therefore 



Idet^^l = n{f4) 2 cos 2 ^ + 4sin 2 ^), (360) 



r=0 



2 AT— 1 _ AT-1 / ^ 



*njdet^| = /3 2 n{l+(^) 2 } • 



(362) 



4sinh 2 -. (363) 



Since the two solutions of (375) have always opposite signs, the Maslov 
index of this part is 0. 

B.3 Inverse hyperbolic type 

In this case, the normal form of V is 
cos 2-7rf — sin 2nt 



, . sin 2irt cos 27rf . 
^ (*) = { / ./5(2*-D n 7 \ ( 364 ) 



(0<*<l/2), 

/ _ e /3(2t-l) o " \ 

( - e -«*-i)J (V2<*<1) 



We calculate the trace of the quantum time evolution operator corre- 
sponding to this V: 

Z = TrM(V(i = l)), (365) 
= Tr {M (V(l, 1/2)) M (V(l/2, 0))} . (366) 
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Here, V(ti,t 2 ) denotes the classical propagator from t = t\ to t = t 2 - 
Let us evaluate this trace using the coherent states \z) = \q + ip). 

Z = ( Zl \M (V(l, 1/2)) \z 2 )(z 2 \M(V(l/2,l))\ Zl l367) 

J IT 7T 

Since (z 2 \M(V (1/2, l))\ Zl ) = (-z 2 \z{)e-™/ 2 , 

/r/z 
-<*|M(V(l,l/2))|-*>. (368) 

Therefore to calculate the path integral of the inverse hyperbolic type 
is equivalent to calculate the path integral of the hyperbolic type with 
anti-periodic boundary condition. The eigenvalue equation is 

Qj+i - Qj + 2jf(.Qj+i + Qj) = cPj+i. ( 369 ) 
Pj-Pj+i + ^iPi+Pj+i) = «& (l<j<N-l). (370) 



<7i + Qn + 2tf(Qi - In) = epi, (371) 
Pn +Pi + TjjjiPN -Pi) = £Qn (372) 
We assume the solutions of eigenvalue equations as 

:;)=(:M iS ^) <'-m.~."-d. 



. ir(2r+l) a , - tt(2j-+1) 

1-e" 1 ^-^ + ^.(i + e -*-h^) 



Then the eigenvalue equation becomes 

( 

■ 7r(2r+l) a , . 7r(2r+l) „ 

\ 1 - e'^rr^ + ^(1 + e-h*- 1 ) 

(374) 

2 f I3 2 2 7r(2r + l) „ . 2 7r(2r + l)\ n 
Therefore 

= n{(!)- 2 ^ + W^}. (are, 
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lim \detp N \ 

N— »oo 



r=0 

= 4 cosh 



1 4 

2 







ir(2r + 1) 



(378) 
(379) 



Here we used Q415| ) and (409). Since the inverse hyperbolic type can 
be deformed continuously to elliptic type (f3 — >0), the Maslov index of 
this part is the same as the elliptic type. Therefore Maslov index of 
this part is 1. 

B.4 Loxodromic type 

The normal Hamiltonian of this type is 

H = a(piq 2 —p2qx) - Pijpiqx + P2Q2)- (380) 
The eigenvalue equation is 





- "+ 




"2iV (g2j " 




= £Pi(381) 


12, j 


- 92j-l 4 








= ^382) 


Phi 


" P2,j+1 + 




a . 


l"P2,j+l) 


= a?i(333) 


P2,j 


- P2,j+1 - 


a . . 


/3 , 
f 2N {P ^' 


l"P2,j+l) 


= ^384) 



If we assume the forms of eigenvectors of Tf> N as 

(r = 0,1, ...,#-!), (385) 



9ij 



f PI \ 

P2 

V 92 J 



2nr 

cxp 1 l W J 



the eigenvalue equation become 



-4,. 



1 — e 1 « 



A r 
At 



#(l + e-^) 







( Pi \ 




= e 


P2 


9i 




9i 


V 92 J 




V 92 ) 



(386) 



2_fl + e -» 

27V V- 1 - ^ e 

o-<^F + JL(\ 



- e 1 jv ) 



det(e 2 / - A r A+) = 0. 



(387) 
(388) 
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The explicit form of A r A\ is 

A A> = ( T cos2 lT+ 4sin2 TNr 4ifsin^cos^ \ 
rr V -4*#sin^cos^ ^cos^+4sin 2 ^ J ' 

(389) 

Therefore the eigenvalue equation becomes 

e4 - 2 (^2^ cos2 77 + 4 sin2 77) f2 + det ^4 = 0, (390) 

det = (4 sm 2 - + cos 2 -j + ^f- cos 4 - (391) 

Hence 

AT-l 

Idet^^l = JJ detAr^, (392) 



r=0 



iV 2 



r=l 
JV-1 



™ 2 nUl + ^ot 2 ^ + ^cot^ . (394) 



Jimjdet^l = (a 2 +(3 2 ) 2 l[ 1 



(27rr) 2 

2 \ 4 00 



(2a/3) 2 
(2?rr) 



(395) 



where 



7 
5 2 



q 2 -/3 2 

(2tt) 2 

2a/3 



Thus we obtain the result 
\det$ N \ 



(27r)2- 
4 (cosh /3 — cos a) 2 , 



(397) 
(398) 

(399) 



using the formulas (407) and (411). 

Since the eigenvalue equation (388) always has two positive and 
two negative eigenvalues, the Maslov index of this part is 0. 
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B.5 Parabolic type 

The Hamiltonian is 

and the eigenvalue equation is 



TT- 7 „ 2 



(400) 



1 





- e » 



--3-(2 



1 - e' 1 ' 
+ e l ~ 



0. 



. (401) 
(402) 



e 2 + 1 

AT TV TV 

If r 7^ 0, the eigenvalues of ( |402j ) are not equal to zero, and the deter- 
minant of this part is 



N 

n 



-4 sin 



N 



(-1) 



N-1 N 2 



(403) 



Since the two solutions of (402) have opposite signs, these eigenvalues 
don't contribute to the Maslov index. 
If r = 0, the eigenvalues are 



(404) 



and the Maslov index of this part is ^sgn7. The normalized eigenvector 
corresponding to these eigenvalues are 



Pi 



(405) 



The volume factor of the zero-mode measured by these vectors are 
Vjv = \/~NV, where V is the volume factor by the normal measure in 
the phase space. Therefore the contribution from this part is 



dX 



N 



(2TTh) 



N 



exp 



4,x T N U)x 



2h 



■N 



v N 



V 



^2-nUN-i V27r?u7 



(406) 
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Useful formulas 

oo , 

n (•-.. 

n=l v 

°° / T 2 



x 2 \ sin 7rx 



sinh: 



,,2 



a; 2 \ cosh 2; — cos a 



n i 1 

r 

00 j 



(a + 2n7r) 2 / 1 — cos a 



cosh 2/ — cos 2i? 



n=l 



(n 2 -a 2 )V V^T^ 1 2 si 



sin 7ra 



i? = Re 7r\/ a 2 + zx 2 
I = Im 7r\/ a 2 + ix 2 



N-l 

n 

r=l 



sin 



7rr TV 



11 sm 



r=0 



2iV 
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